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ABSTRACT 


Very recently, molecular rotational transitions observed in the photon- 
dominated region of the Horsehead nebula have been attributed to /-C 3 HV In 
an effort to corroborate this finding, we employed state-of-the art and proven 
high-accuracy quantum chemical techniques to compute spectroscopic constants 
for this cation and its isotopologues. Even though the B rotational constant from 
the fit of the observed spectrum and our computations agree to within 20 MHz, 
a typical level of accuracy, the D rotational constant differs by more than 40%, 
while the H rotational constant differs by three orders of magnitude. With the 
likely errors in the rotational transition energies resulting from this difference in 
D on the order of 1 MHz for the lowest observed transition (J = 4 — > 3) and 
growing as J increases, the assignment of the observed rotational lines from the 
Horsehead nebula to Z-C 3 H+ is questionable. 

Keywords: astro chemistry — ISM: individual objects: Horsehead nebula — ISM: lines 
and bands — ISM: molecules — molecular data — radio lines: ISM 


1. Introduction 

Within the past months, /-CaH" 1 " has been attributed as the source of several rotational 
lines observed in the Horsehead nebula photon-dominated region (PDR) (Pety et al. 2012). 
In order for the observed spectral features to be linked to a single molecular carrier, the 
molecule must be linear with a ground state and a lone (R-type) rotational constant close 
to 11.24 GHz. However, only standard quantum chemical computations of Z-CsH + have 
given qualitative corroborating evidence for the rotational constants, and no quant it at ive, 
high-accuracy reference data, either from computation or laboratory experiment, exists for 
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Z-C 3 H+. The corresponding C 3 H radical has long been known to exist in the interstellar 
medium (ISM) (Thaddeus et al. 1985) and was confirmed with concurrently generated 
laboratory reference data (Gottlieb et al. 1985). Although laboratory studies of known 
interstellar molecules in simulated interstellar environments have often produced what is 
believed to be C 3 H + (Bohme et al. 1983; Schwell et al. 2012) with theory providing further 
insights (Radorn et al. 1976; Ikuta 1996; Wang et al. 2007), neither high-accuracy rotational 
nor vibrational spectroscopic data has yet been generated for Z-C 3 H + . 

Recently, we have employed state-of-the art ab initio quartic force fields (QFFs) 
and perturbation theory at second order (Papousek & Aliev 1982) to compute rotational 
constants that are within 20 MHz or better of experiment (Fortenberry et al. 2011a, b, 
2012b, a) for similar tetra-atomic systems. The use of QFFs has also been applied to 
other molecules of interstellar relevance in order to compute accurately their spectroscopic 
constants in addition to their fundamental vibrational frequencies and even some overtones 
and combination bands (Huang & Lee 2008; Lee et al. 2009; Huang & Lee 2009; Inostroza 
et al. 2011; Huang et al. 2011; Huang & Lee 2011; Fortenberry et al. 2012a, 2013b). 
Since high-accuracy reference data is currently not available for Z-C 3 H + , the methods 
utilized previously are employed here in order to assist in confirmation for the detection 
of this cation in the Horsehead nebula. Additionally, we will also supply spectroscopic 
constants and vibrational frequencies for the deuterated form and for each of the three 13 C 
singly-substituted isotopologues. 


2. Computational Details 

Utilizing the MOLPRO 2006.1 suite of quantum chemical programs (Werner et al. 
2006), spin-restricted Hartree-Fock (RHF) (Schemer et al. 1987) coupled cluster theory 
(Lee & Scuseria 1995; Shavitt & Bartlett 2009; Crawford & Schaefer 2000) at the singles, 
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doubles, and perturbative triples level [CCSD(T)] (Raghavachari et al. 1989) is employed 
to compute the reference geometry. This initial geometry is based on a composite approach 
to give the most accurate bond lengths in this linear molecule. The CCSD(T) bond lengths 
(denoted as R) computed with the cc-pV5Z basis set (Dunning 1989; Peterson & Dunning 
1995; Dunning et al. 2001) are corrected for effects from the smaller cc-pVQZ basis set as 
well as core correlation and scalar relativistic effects (Douglas & Kroll 1974) through the 
following approach: 


R — RbZ + (R5 Z — RqZ ) + ( Rnel.+Core ~ RqZ )• (1) 

The Rpiei.+Core CCSD(T) geometry optimization makes use of the standard cc-pVQZ-DK 
basis set, which can account for relativistic effects, further augmented to include core 
correlating s, p, d , and / functions from the standard cc-pCVQZ basis set (Dunning 1989). 
The last term in Eq. 1 serves to isolate the scalar relativistic and core correlating effects so 
that they may then enhance the standard cc-pV5Z bond lengths. 

From this reference geometry, the quartic force field (QFF) is created from 
displacements of 0.005 A and 0.005 radians (Huang et al. 2011), respectively for the three 
bond lengths and two doubly-degenerate linear bends, and contains 569 points. At each 
point CCSD(T)/cc-pVTZ, cc-pVQZ, and cc-pV5Z energies are computed and extrapolated 
to the complete basis set (CBS) limit via a three-point formula (Martin & Lee 1996). 
The CBS energy is then corrected for scalar relativity (Douglas & Kroll 1974) and core 
correlation effects from the Martin-Taylor basis set (Martin & Taylor 1994). This results 
in the CcCR QFF defined previously (Fortenberry et al. 2011a, 2013a). The QFF is then 
fit with a sum of residual squares of 5.7 x 10~ 1( a.u. 2 to a Taylor series expansion of the 
potential and is of the form: 

1=2 'y 1 -FijAjAj + — 'y " KjfcjAjAj Afc + — N ^ FAjzAjAj A^A/ (2) 

ij ijk ijkl 
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where F t j are force constants and Aj are the displacements. This high-symmetry C^, 
system is most easily represented in the simple-internal coordinate system defined here 
(with atom numbering from Fig. 1) as: 


Si(£ + ) 

= Cl -H 

(3) 

s 2 (e + ) 

= Ci - c 2 

(4) 

s 3 (£ + ) 

= C 2 - c 3 

(5) 

5 4 (n„) 

= LIN1(H - Ci - C 2 -y) 

(6) 

s 5 ( n yz ) 

= LIN1(H - Ci - C 2 - x) 

(7) 

s 6 {n X z) 

= LINl(Ci — C 2 — C 3 — $ 

(8) 

SriM 

= LINl(Ci - C 2 - C 3 - x) 

(9) 


These coordinate subscripts correspond to the force constants given in Table 1. The x and 
y quantities are necessary to define the linear bending angles (S 4 — S 7 ) and correspond to 
a direction perpendicular to the bending plane so that degenerate linear modes may be 
uniquely specified, since x -L y. Second-order vibrational perturbation theory (VPT2) (Mills 
1972; Watson 1977) via the SPECTRO program (Gaw et al. 1996) produces the desired 
spectroscopic constants and vibrational frequencies. Since the QFF is computed within the 
Born-Oppenhcimer approximation, only the SPECTRO anharmonic analyses differ between 
the five total isotopologues. 


3. Discussion 

The /-C 3 H + and isotopologue bond lengths, rotational constants, quartic and sextic 
centrifugal distortion constants, and center-of-mass dipole moments are given in Table 2, 
and the fundamental vibrational frequencies are listed in Table 3. The relevant data are 
in fairly good agreement with the full-dimensional PES computed previously given the 


5 


much lower level of theory used (Wang et al. 2007). The inclusion of the scalar relativity 
correction is necessary since, for example, B 0 for standard Z-C 3 H + computed with only 
the CBS-energy-including QFF is 11 209.78 MHz, while the CR QFF (which additionally 
includes the relativistic terms, see Fortenberry et al. (2011a)) raises this value to 11 213.99 
MHz. Similarly but more importantly, core correlation is also necessary as the rotational 
constant for the CcC QFF, which is the CBS plus core correlation QFF, puts Bo at 11 
258.41 MHz. The total CcCR QFF Bq, listed in Table 2, is thus 11 262.68 MHz. This 
value differs by 17.73 MHz from that (11 244.947 4 MHz) obtained in the fitting by 
Pety et al. (2012) to assign the observed lines in the Horsehead nebula to /-C 3 HQ This 
represents similar accuracy as that noted previously for the cis- and frans-HOCO radicals 
(Fortenberry et al. 2011a, b), frans-HOCO + (Fortenberry et al. 2012b), and frcms-HOCS + 
(Fortenberry et al. 2012a), as well as many other examples. 

However, the correspondence between theory and experiment for the H-type constant 
is not as good. The CcCR D e value is 4.248 kHz while the D value obtained by Pety 
et al. (2012) in their second-order fit is 7.652 kHz and 7.766 kHz for their third-order fit. 
This represents a difference of 3.404 kHz or 44.5% between theory and experiment for the 
second-order fit D and is larger still for the third-order fit. Though this is a comparison of 
D e to D 0 , i.e., vibrational averaging is not included in the theoretical value, the percentage 
change between the computed Bo and B e values for this same molecule is only 0.77%. 
Further, various QFFs have been formulated to examine this system (33 in all) ranging from 
those that are state-of-the art and have been proven to yield highly-accurate spectroscopic 
constants including D , such as the given CcCR QFF, to those that use more modest levels 
of theory, such as CCSD(T)/cc-pVTZ, and the range in D e values for these QFFs is only 
0.053 kHz, which conclusively shows that the ab initio D e value is converged to an extent 
that rules out a value near 7.6-7 . 8 kHz. 
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Additionally, other linear systems have been studied using similar theoretical 
approaches. For instance, the percentage error between the computed D e (84.135 kHz) and 
the experimental D [87.22 kHz from Winnewisser et al. (1971)] is 3.5% for HCN computed 
with a CCSD(T)/AN01 QFF (Lee et al. 1993). The D/D e error for CCH~ is 2.1% as 
computed with a CcCRE QFF (with corrections for higher-order electron correlation effects: 
“E”) in the same second order perturbational approach utilized here, where D e is 94.3 
kHz (Huang & Lee 2009), and the experimental D is 96.97 kHz (Briinken et al. 2007). 
Furthermore, the error for acetylene is only 2.27% between a D e of 47.673 kHz (Martin 
et al. 1998) and a D of 48.780 kHz (Kabbadj et al. 1991). Perhaps a more relevant example 
for the current case is given for C3N - where Do has been experimentally determined to 
be 0.68 578 kHZ and a CCSD(T)/aug-cc-pVQZ QFF yields a D e of 0.627 kHz (Thaddeus 
et al. 2008; McCarthy & Thaddeus 2008). The 8.57% error is still much less than the 
44.5% discrepancy for CsH + and is probably somewhat larger then the other examples 
because of using lower levels of theory to compute the QFF, and because it is well known 
that small molecular anions are harder to describe accurately. Hence, the known errors 
for these systems are an order of magnitude smaller than that which is present for D e and 
D in Z-C 3 H+. As a result, we cannot attribute the substantial disagreement between the 
computed D e for Z-C 3 H + and the D value obtained by Pety et al. (2012) merely to errors 
in the theoretical approach. Thus, it is doubtful that the series of lines observed in the 
Horsehead nebula actually correspond to Z-C3HA 

If the computed D e value closely represents the actual D for /-C3H + and assuming 
the Bq value reported by Pety et al. (2012), the J = 4 — y 3 transition will differ from that 
reported by Pety et al. (2012) by 0.9 MHz. The J = 5 — » 4 transition will differ by 1.7 
MHz, J = 6 — y 5 by about 3.0 MHz, and these differences will continue to increase up to 
23.4 MHz as J increases to the J = 12 — y 11 transition. For high-resolution rotational 
spectroscopy, these large discrepancies are well outside of the precision present in most 
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spectrometers utilized in the laboratory and for observation of the ISM, and well outside 
the errors expected from state-of-the art ab initio calculations. For example, Brites & 
Jutier (2012) computed the rotational transitions of HCNH + with a six-dimensional PES 
from a CCSD(T)-F12A/cc-pCVTZ-F12 QFF for J = 1 — > 0 up to J = 10 — > 9 with errors 
ranging from 0.6 MHz to 4.84 MHz compared to experimental work (Amano et al. 2006). 
This comparison of theory and experiment again strongly supports our assertion that the 
assignment of the observations in the Horsehead Nebula PDR to Z-CsH + are incorrect. 

Similarly, the computed H e value for Z-CsH + is 0.375 mHz, while the third-order fit 
H obtained by Pety et al. (2012) is 0.56 Hz or 560 mHz. This represents three orders of 
magnitude difference between the two values. There are not as many comparisons in the 
literature between theoretical H e values and experimental H values, but for acetylene, H e 
has been computed to be 0.0380 Hz (Martin et al. 1998), while H has been experimentally 
determined to be 0.0480 Hz (Kabbadj et al. 1991), an error of 20.8%. Even if the computed 
H e for Z-C 3 H + is in error from the proper physical value by as much as 40%, this is only a 
minute step towards the value necessary to fit the lines observed. Admittedly, comparison 
of computed H e values to experimental H constants has not been as common making 
the error range for the computed value less certain, but the substantially large difference 
between values here should be well-beyond the potential accuracy range for the computed 
H e . This difference in the H constant casts further doubt on /-CsH + as the carrier of the 
observed transitions in the Horsehead nebula PDR. 

Since the energy points necessary to define the QFF are computed within the 
Born-Oppenhcimer approximation, the same QFF can also be used to produce VPT2 
spectroscopic constants and vibrational frequencies for the isotopologues, as well. Thus, the 
spectroscopic constants of C 3 D + , 13 CCCH + , C 13 CCH + , and CC 13 CH + are listed in Table 
2, and the fundamental vibrational frequencies are given in Table 3. These reference data 


should aid in the analysis of further laboratory studies of interstellar reactions that create 
Z-C 3 H+. 


4. Conclusions 

Even though the CcCR VPT2 computed B 0 rotational constant for Z-C 3 H + differs 
by less than 20 MHz from that obtained by Pety et al. (2012) to assign lines in the 
Horsehead nebula PDR to this cation, the more substantial differences between the D and 
H constants strongly questions the assignment of /-C 3 H + to the observed lines. These large 
discrepancies, especially for the D constant, will alter the rotational spectrum of X 1 E + 
C 3 H+ by nearly 1 MHz for the lowest observed rotational transition, J — 4 — > 3, and by 
more than 23 MHz for the highest, J = 12 — )■ 11. Hence, the astronomical lines observed 
cannot be conclusively linked to this cation as the carrier. 

It is hoped that the reference data provided here should assist in laboratory studies 
that can conclusively determine the rotational spectrum of /-C 3 H + . Additionally, the data 
are extended to the deuterated and 13 C singly-substituted isotopologues as well. C 3 H + has 
long been hypothesized to exist in the ISM and is of significance in the carbon chemistry 
of space (Turner et al. 2000; Wakelam et al. 2010). There is little doubt that it should be 
present in ISM, but the attribution of the observed lines by Pety et al. (2012) to Z-C 3 H + is 
questionable. 
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Tabic 1: The simple-internal CsH + CcCR QFF Force Constants (in rndyn/A n -rad m ). 


Fn 

5.810 273 

F441 

-0.1688 

F2222 

405.81 

F5532 

-0.16 

F7511 

0.28 

F21 

-0.214 189 

F442 

-0.6565 

F3111 

0.05 

F5533 

0.22 

F7521 

0.94 

F22 

13.840 846 

F443 

-0.1210 

F3211 

0.22 

t? a 

^ 5544 

0.41 

F7522 

-0.11 

F31 

0.075 121 

F551 

-0.1688 

F3221 

0.22 

F5555 

1.22 

F7531 

0.01 

F32 

0.116 818 

F552 

-0.6565 

F3222 

2.09 

F 6411 

0.28 

F7532 

0.36 

F33 

8.520 119 

F553 

-0.1210 

F3311 

0.15 

F 6421 

0.94 

F7533 

-0.35 

F44 

0.330 792 

F 641 

-0.0004 

F3321 

0.30 

F 6422 

-0.11 

F 7544 ° 

-0.04 

F55 

0.330 792 

F 642 

0.0412 

F3322 

-1.51 

F 6431 

0.01 

F 7555 

-0.11 

F 64 

0.004 472 

F 643 

0.1390 

F3331 

-0.05 

F 6432 

0.36 

F 7654 a 

-0.05 

F66 

0.034 854 

F 661 

-0.0212 

F3332 

-0.98 

F 6433 

-0.35 

t? a 

^ 7665 

-0.03 

F75 

0.004 472 

F 662 

-0.2630 

F3333 

245.44 

F 6444 

-0.11 

F7711 

-0.08 

F77 

0.034 854 

F 663 

0.0347 

F4411 

-0.08 

T 7 ' a 

^ 6554 

-0.04 

F7721 

0.10 

F111 

-33.4159 

F75I 

-0.0004 

F4421 

0.45 

F6611 

-0.08 

F7722 

0.18 

F211 

0.1250 

F752 

0.0412 

F4422 

-0.18 

F 6621 

0.10 

F7731 

0.11 

F221 

0.0463 

F753 

0.1390 

F4431 

-0.22 

F 6622 

0.18 

F7732 

-0.38 

F222 

-82.4686 

F77I 

-0.0212 

F4432 

-0.16 

F 6631 

0.11 

F7733 

-0.19 

F311 

-0.0121 

F772 

-0.2630 

F4433 

0.22 

F 6632 

-0.38 

F7744 

0.00 

F321 

-0.0241 

F773 

0.0347 

F4444 

1.22 

F 6633 

-0.19 

F7755 

-0.10 

F322 

-0.7100 

Full 

172.62 

F5511 

-0.08 

F 6644 

-0.10 

F 7764 ° 

-0.03 

F331 

-0.1360 

F2III 

0.14 

F5521 

0.45 

F 6655 

0.00 

T? a 

7766 

0.02 

F332 

-0.4027 

F22II 

-0.69 

F5522 

-0.18 

F 6664 

-0.10 

F7775 

-0.10 

F333 

-49.4860 

F222I 

0.98 

F5531 

-0.22 

F 6666 

0.06 

F7777 

0.06 


“These necessary force constants are not symmetry-unique but are defined from relationships of 
other force constants. See Ref. 60 of Martin et al. (1998). 
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Table 2: The CcCR QFF Zero-Point (R a vibrationally-averaged) and Equilibrium Struc- 
tures, Rotational Constants, CCSD(T)/aug-cc-pV5Z Dipole Moment, Vibration-Rotation 
Interaction Constants, and Quartic and Sextic Distortion Constants of linear X 1 E + C 3 H + , 
C 3 D+, 13 CCCH+, C 13 CCH+, and CC 13 CH+. 



c 3 H+ 

C 3 D+ 

1J CCCH+ 

C id CCH+ 

CC 1;i CH+ 

ro(Ci-H) 

1.069 951 A 

1.072 922 A 

1.067 560 A 

1.070 095 A 

1.069 996 A 

r 0 (Ci-C 2 ) 

1.230 766 A 

1.231 576 A 

1.230 543 A 

1.230 961 A 

1.230 903 A 

ro(C 2 -C 3 ) 

1.332 581 A 

1.331 723 A 

1.332 879 A 

1.332 800 A 

1.332 416 A 

Bo 

11 262.68 MHz 

10 115.99 MHz 

10 824.96 MHz 

11 258.19 MHz 

10 919.88 MHz 

a B 1 

33.5 MHz 

49.0 MHz 

31.5 MHz 

33.2 MHz 

30.6 MHz 

a B 2 

74.8 MHz 

55.0 MHz 

71.5 MHz 

72.2 MHz 

73.4 MHz 

a B 3 

37.1 MHz 

30.1 MHz 

35.4 MHz 

37.2 MHz 

35.1 MHz 

a B 4 

-2.6 MHz 

-10.2 MHz 

-2.3 MHz 

-2.5 MHz 

-1.7 MHz 

ct B 5 

-157.3 MHz 

-133.2 MHz 

-151.2 MHz 

-157.7 MHz 

152.7 MHz 

1- a 
> aaaa 

-16.992 kHz 

-12.991 kHz 

-15.803 kHz 

-16.988 kHz 

-15.966 kHz 

r e (Ci-H) & 

1.078 961 A 

- 

- 

- 

- 

r e (Ci-C 2 ) 

1.235 360 A 

- 

- 

- 

- 

r e (C 2 — C 3 ) 

1.339 841 A 

- 

- 

- 

- 

B e 

11 175.51 MHz 

10 039.77 MHz 

10 740.68 MHz 

11 175.30 MHz 

10 834.97 MHz 

D e 

4.248 kHz 

3.248 kHz 

3.951 kHz 

4.247 kHz 

3.992 kHz 

H e c 

0.375 mHz 

0.315 mHz 

0.315 mHz 

0.375 mHz 

0.361 mHz 

d 

H'Z 

1.20 D 

- 

- 

- 

- 


“Since this is a linear molecule, r a aaa = Tbbbb = r aa {,5 while all other quartic centrifugal distortion constants are 0.0. 
& The use of the Born-Oppenheimer approximation necessitates that the equilibrium geometry be the same for each 

isotopologue. 

c H e = $aaa = &bbbi and all other sextic centrifugal distortion constants are 0.0. 
d The C3H+ coordinates (in A with the center-of-mass at the origin) used to generate the Born-Oppenheimer dipole moment 
component are: H, 0.000000, 0.000000, -2.285175; Ci, 0.000000, 0.000000, -1.206214; C 2 , 0.000000, 0.000000, 0.029147; C 3 , 

0.000000, 0.000000, 1.368989. 


Table 3: The VPT2 CcCR QFF harmonic and anharmonic fundamental vibrational frequen- 
cies (in cnr 1 ) for C 3 H+, C 3 P+, 13 CCCH+, C 13 CCH+, and CC 13 CH+. 


Mode 

Description 

C 3 H+ 

Harm. Anharm. 

C 3 D+ 

Harm. Anharm. 

la CCCH+ 
Harm. Anharm. 

C ia CCH+ 
Harm. Anharm. 

CC l;i CH+ 
Harm. Anharm. 

IZl(o-) 

Ci— H stretch 

3309.7 

3167.8* 

2580.0 

2502.2 

3309.7 

3166.7* 

3308.2 

3163.7* 

3293.1 

3152.0* 


C 2 — C3 stretch 

2142.7 

2096.3* 

2015.5 

1972.2 

2136.4 

2090.4* 

2088.4 

2044.9* 

2126.1 

2080.0* 

v 3(cr) 

Ci— C 2 stretch 

1189.3 

1194.1* 

1162.8 

1162.2 

1161.3 

1167.2* 

1188.5 

1191.2* 

1172.7 

1178.5* 

l/ 4 (ff) 

H-Ci-C 2 bend 

805.8 

782.3 

641.7 

625.8* 

805.8 

782.2 

804.3 

781.0 

799.0 

776.0 

t'Sp) 

Ci-C 2 -C 3 bend 

124.0 

114.2 

117.8 

110.6* 

123.2 

113.4 

121.0 

111.2 

123.7 

114.1 


* Denotes states in Fermi resonance. 
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Fig. 1. — The CcCR equilibrium geometry of X 1 £ + CsH + . 
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